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COMPUTER  PROGRAMS  FOR  HAYSTACK  SERVO  TESTING 


1.  INTRODUCTION 

It  is  intended  that  the  Haystack  antenna  be  primarily  directed  by  the  490  digital 
computer.  When  the  digital  control  system  is  first  connected  to  the  antenna  and 
actual  motion  is  attempted,  it  will  be  important  to  study  the  response  characteristics 
of  the  antenna  servo  system.  It  will  also  be  desirable  to  periodically  check  this  re¬ 
sponse  over  the  useful  life  of  the  facility.  Although  the  system  response  can  be 
checked  by  analog  techniques,  and  although  the  normal  ^'Pointing”  program  also  pro¬ 
vides  some  checking,  it  was  felt  that  a  special  effort  was  justified.  Specifically,  it 
appeared  desirable  to  permit  digital  generation  of  special  test  functions  and  a  detailed 
digital  analysis  of  the  resulting  response.  For  this  purpose,  a  computer  program  has 
been  written  for  the  Univac  490  which  can  drive  the  antenna  and  which  will  provide  on¬ 
line  printouts  in  real  time  of  the  input  to  the  servo  and  the  response  of  the  servo.  For 
more  detailed  checking,  the  program  provides  magnetic  tape  outputs  which  can  be 
further  processed  and  plotted  using  the  Laboratory  IBM  7094  in  Lexington. 

Probably  the  most  informative  viewpoint  to  take  in  attempting  to  comprehend  all 
the  ramifications  of  this  system  of  programs  is  to  consider  it  as  a  completely  flexible 
and  automatic  data-gathering  system.  It  is  not  the  function  of  the  system  to  provide 
analyses  of  complex  aspects  of  antenna  steering  kinematics.  The  function  of  the 
system  is  to  provide  a  perfectly  general,  systematic,  and  uniform  procedure  for 
automatically  obtaining  static  or  dynamic  data  about  antenna  behavior  under  a  variety 
of  conditions.  The  data  obtained  is  in  printed,  magnetically  recorded  (digital),  or 
graphed  form  and  is,  therefore,  suitable  for  immediate  analysis  by  manual  or  mechani¬ 
cal  computation  procedures.  In  particular,  the  magneitc  tape  records  are  acceptable 
to  the  IBM  7094  computer  for  further  computation. 

This  testing  is  performed  by  entering  through  the  on-line  console  specifications 
for  independent  motions  to  be  executed  by  the  azimuth  and  elevation  drive  systems. 

One  is  free  to  specify  dynamic  tests  with  sinusoidal  oscillations,  step  and  impulse 
functions,  and  up  to  fifth  order  polynomial  trajectories.  A  static  positioning  test  is 
also  available.  The  antenna  servo  response  to  these  signals  is  automatically  recorded 
on  magnetic  tape  and  simultaneously  printed  on  the  on-line  printer.  The  magnetic 
tape  can  be  further  processed  by  the  IBM  7094  system  to  produce  graphs  of  the  re¬ 
sponse  (via  the  Calcomp  plotter)  and  BCD  tapes  suitable  for  printing  or  input  to 
FORTRAN  programs. 
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The  remainder  of  this  report  discusses  the  testing  package  in  increasingly  greater 
detail.  First,  typical  examples  of  the  data  generated  are  exhibited.  Detailed  instructions 
for  system  usage  follow,  terminating  in  a  detailed  documentation  of  the  mechanical 
logic  of  the  computer  programs.  Examples  and  discussions  tend  to  be  in  terms  of  the 
particular  aspect  of  the  problem  with  which  our  group  is  associated;  however,  this 
aspect  is  extremely  general  itself. 

II.  SYSTEM  DESCRIPTION 

In  using  the  system,  one  first  performs  certain  routine  operations  which  bring  the 
program  into  the  Univac  490  computer.  When  the  program  is  started,  it  makes  a  few 
basic  inquiries  and  then  leaves  the  operator  free  to  specify  his  requirements.  The 
specifications  of  a  run  permit  complete  freedom  in  selecting  starting  positions  (biases) 
on  both  axes.  Nominally,  the  experiment  begins  at  a  specified  time  and  lasts  for  a 
specified  duration,  with  motions  commencing  when  the  antenna  has  stabilized  at  the 
specified  starting  point.  One  always  has  full  and  independent  control  over  the  amount 
of  data  recorded  on  magnetic  tape  and  on  the  on-line  printer.  The  mathematical 
functions  which  define  the  drive  functions  (see  below)  can  be  independently  specified  on 
both  axes.  An  experiment  can  be  interrupted  at  any  time  without  losing  any  data  previ¬ 
ously  obtained. 

Five  driving  functions  are  available  on  each  axis.  AZNOT  and  ELNOT  simply  hold 
the  antenna  fixed  at  a  specified  azimuth  or  elevation  bias.  AZSINE  and  ELSINE  produce 
a  sinusoidal  oscillation  about  the  bias  angle  with  independently  specified  amplitude  (in 
degrees),  and  frequency  (in  cycles  per  second).  AZPOLY  and  ELPOLY  cause  the 
antenna  to  track  any  trajectory  that  can  be  specified  parametrically  in  time  by  an 
equation  of  the  fifth  order,  or  less. 

AZIMUTH  =  a  +  a  ,t  +  a„t^  +  a.t^  +  a  +  a.t^  +  BIAS 
o  I  2  3  4  5 

Positive  or  negative  values  may  be  entered  for  the  coefficients. 

To  test  transient  response,  two  types  of  functions  are  available.  AZSTEP  and 
ELSTEP  provide  a  positive  or  negative-going  staircase  function  of  time,  starting  at 
the  specified  bias  =  STEP(A,  B,  C)  +  BIAS: 
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angle 
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1 

1 

B 

i 

1 

^  ^ 

BIAS 

^  w 

1 

^ -  A - ► 

- :  9 

time 

The  parameter  B  may  be  negative,  thus  providing  the  negative -going  function.  The 
other  test  available  is  provided  by  AZIMPULSE  and  ELIMPULSE  which  produce  a 
positive  or  negative  directed  series  of  impulse,  also  beginning  at  the  specified  bias: 
IMPULSE  (A,B,C,D)+  BIAS: 


angle 


Again,  the  parameter  B  may  be  positive  or  negative. 

Once  an  experiment  has  been  specified  by  the  operator,  the  Univac  490  computer 
program  generates  instructions  for  the  antenna  servo  system  at  the  standard  rate  of 
250  times  per  second.  The  data  is  transferred  to  the  servo  system  in  exactly  the  same 
fashion  as  the  data  is  handled  in  the  standard  "Pointing"  program. 


3 


In  real  time,  two  forms  of  output  are  produced;  hard  copy  on  the  on-line,  higli- 
speed  printer  and  IBM  format  magnetic  tapes  for  further  processing  at  the  7094.  It  is 
anticipated  that  the  line  printer  output  will  be  used  for  more  simple  testing,  but  that 
the  detailed  analysis  may  be  useful  for  trouble  analyses  or  for  graphical  presentation 
of  results.  Due  to  the  speed  limitations  on  the  on-line  printer,  it  is  not  possible  to 
obtain  by  that  method  the  full  rate  of  250-point-per-second  data;  the  maximum  output 
rate  of  the  line  printer  in  this  program  is  10  lines  per  second  (or  a  sampling  of  the 
servo  input  and  output  at  1/25  of  the  actual  data  transfer  rate).  The  magnetic  tape 
recordings,  of  course,  provide  complete  data  on  every  point  in  the  full  rate. 

The  magnetic  tapes  generated  by  these  tests  contain  a  record  of  the  signal  applied 
to  the  servos  and  the  position  of  the  axes  when  the  signal  was  applied,  along  with  the 
time  of  application.  A  special  program  on  the  7094  will  convert  these  tapes  to  a  BCD 
format  suitable  for  printing  or  input  to  another  7094  program.  Another  7094  program 
produces  graphs  of  selected  material  directly  from  these  tapes  with  completely  labeled 
and  calibrated  axes.  These  graphs  are  drawn  by  the  Calcomp  plotter  via  the  IBM  1401. 

III.  OPERATING  INSTRUCTIONS  FOR  HAYSTACK  TESTING  PACKAGE 

The  testing  facility  for  the  Haystack  antenna  servo  system  consists  of  three  pro¬ 
grams: 

Program  I  to  run  on  the  Univac  490  computer  which  supplies  commands  to  the 

servos  and  records  their  responses. 

Program  II  to  run  on  the  IBM  7094  to  convert  and  print  the  output  of  Program  I. 

Program  III  to  run  on  the  IBM  7094  to  plot  the  output  of  Program  I  on  the  Calcomp 

plotter. 

A.  Program  I  -  Univac  490  Testing  Program 

The  following  is  an  exact  explanation  of  one  way  in  which  the  antenna  testing  pro¬ 
gram  (ANTENATEST)  may  be  operated  so  as  to  take  full  advantage  of  all  its  features. 
This  does  not  represent  the  full  generality  of  the  system  since  particular  operations 
are  here  prespecified  when  a  number  of  different  but  effectively  equivalent  alternatives 
exist.  Accordingly,  this  should  prove  particularly  useful  to  those  persons  who  are 
relatively  unacquainted  with  the  testing  system  and/or  the  Univac  490. 

1.  Mount  the  following  magnetic  tapes  on  drives  as  indicated: 
a.  (Logical  0)  SPURT  III  system  tape 
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b.  (Logical  1)  Any  blank  tape 

c.  (Logical  2)  The  ANTENATEST  321  tape 

2.  Press  the  ’^CHANNEL  CLEAR”  buttons  on  the  tape  adapter  cabinet  and 
on  the  printer  control  cabinet. 

3.  On  the  maintenance  panel  of  the  computer,  depress  the  toggle  switches 
marked  ’'OPERATION,  ”  "MASTER  CLEAR,  "  and  "BOOTSTRAP"  in  that  order. 

4.  The  console  printer  will  type  "READY.  "  You  type; 

LD‘~S2l:1  3210  10000  s'^ 

The  character  "o  is  marked  "SPEC"  on  the  keyboard. 

5.  Next,  the  console  printer  will  again  type  "READY.  "  You  type: 

PS  Cl  10000  s  . 

6.  The  console  printer  types: 

THIS  IS  THE  ANTENNA  TESTING  PROGRAM. 

DO  YOU  WISH  TO  HAVE  OUTPUT  ON  THE  ON-LINE  PRINTER. 

IF  SO,  HOW  OFTEN. 

You  may  reply: 

NO  SI 

if  you  wish  to  have  no  on-line  printing,  or 
YES  i: ‘XXX  I  S; 

where  "XXX"  is  a  decimal  integer  of  exactly  three  digits  in  the  range  999  to  025.  If 
you  type  in  a  number  N,  exactly  I/Nth  of  the  data  generated  in  positioning  the  antenna 
will  be  printed  on  the  on-line  printer.  The  output  may  be  obtained  in  octal  or  decimal 
form.  The  decimal  listing  contains  seven  columns.  See  Fig.  2.  From  left  to  right, 
the  first  column  is  time,  in  seconds,  from  the  beginning  of  the  test.  The  second  is  the 
azimuth  position  generated  by  the  computer.  The  third  is  the  current  azimuth  position 
of  the  antenna.  The  fourth  is  the  difference  between  these  figures.  The  fifth,  sixth, 
and  seventh  repeat  these  functions  for  the  elevation  coordinate.  All  angles  are  printed 
in  degrees,  truncated  to  four  decimal  places.  In  octal  mode,  9  columns  are  printed. 
See  Fig.  3.  The  first,  second,  and  fifth  are  identical  to  those  in  decimal  mode.  The 
third  column  is  the  azimuth  command  to  the  antenna  in  octal.  The  fourth  column  is 
the  azimuth  position  of  the  antenna  in  octal.  The  sixth  and  seventh  repeat  the  third 
and  fourth  for  elevation.  The  eighth  and  ninth  are  the  octal  differences  in  azimuth  and 
elevation  between  the  command  and  the  position.  Printing  mode  may  be  selected  in 
mid-run  or  changed  by  typing  the  characters  "O"  for  octal  mode  and  "D"  for  decimal 
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mode.  Printing  may  be  stopped  and  restored  in  mid-run  by  typing  the  characters  "S" 
for  stop  and  "P"  for  print.  None  of  these  characters  will  be  printed  at  the  console. 

7.  The  console  printer  types 

TYPE  TEST  IDENTIFICATION 

You  may  type  a  label  up  to  30  characters  in  length  and  terminated  by  a  (^i  . 

8.  Next,  the  console  printer  types 

TYPE  TEST  PARAMETERS 

You  have  two  options.  You  may  type  ".  .  "  (option  a),  or  you  may  type  a  parameter 

line  (option  b). 

9  a.  If  you  typed  ".  .  ,  "  the  console  will  ask  two  questions.  First,  it 

asks; 


DO  YOU  WISH  TO  CONTINUE 

This  should  be  answered  YES  ^  or  NO  .  Next  it  will  say 
DO  YOU  WISH  TO  REWIND  OUTPUT  TAPE 

Again,  answer  YES  or  NO  .  If  you  answered  YES  to  the  second  question,  the 
tape  on  B  will  be  rewound  with  interlock  and  must  be  changed.  If  you  answered  YES  to 
the  first  question,  you  will  be  back  at  step  6.  If  you  typed  NO,  you  will  be  back  at 
step  5. 

9  b.  If  you  elect  to  enter  data  for  a  new  run,  it  must  be  in  the  form  described 
at  the  end  of  this  discussion.  As  soon  as  you  type  at  the  end  of  the  line,  the  actual 
test  begins.  If  you  enter  an  inadmissible  line,  an  appropriate  remark  will  be  typed  and 
you  will  be  returned  to  step  8. 

10.  After  a  15-second  delay  for  aiming,  the  computer  will  generate  the  speci¬ 
fied  function  and  supply  it  to  the  antenna.  After  specified  delays,  data  will  be  printed 
on-line  and  written  on  magnetic  tape. 

11.  If  you  decide,  in  mid-run,  to  enter  a  new  set  of  parameters  before  the 
full  run  is  complete,  type  .  You  will  be  returned  to  step  7.  All  data  recorded  on 
tape  to  this  point  will  be  usable.  If  you  let  the  run  continue  for  its  specified  duration, 
you  will  still  end  up  at  step  7. 


Note:  It  is  important  to  plan  your  runs  to  take  maximum  advantage  of  identical  constants 
in  printing  rate,  since  the  only  way  to  change  the  rate  is  by  the  jump  from  step  9  a.  to 
step  6,  a  fairly  time  consuming  operation. 
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Format  for  a  Line  of  Parameters 


1.  Type  the  azimuth  bias  in  degrees  as  a  signed  or  unsigned  decimal  number 
(fractional  places  allowed). 

2.  Type  a  comma. 

3.  Repeat  1  and  2  for  elevation  bias. 

4.  Type  the  number  of  seconds  delay  desired  before  recording  output  as  a 
positive  decimal  integer.  Zero  is  not  acceptable.  Type  a  comma. 

5.  Repeat  4  for  the  duration  of  run  desired. 

6.  Type  the  number  of  4  ms  intervals  between  successive  tape  recorded  data 
points  as  a  positive  decimal  integer.  (5  intervals  =  20  ms  between  recorded  measure- 
ments) 

7.  Type  the  name  desired  for  azimuth  drive  function.  It  must  be  AZNOT, 
AZPOLY,  AZSINE,  AZSTEP,  or  AZIMPULSE.  Type  a  comma. 

8.  Repeat  7  for  elevation.  The  names  are  ELNOT,  etc. 

9.  Type  the  arguments,  separated  by  commas,  for  the  azimuth  function. 

10.  Repeat  9  for  the  elevation  function.  At  this  point,  every  field  in  the  line 
should  end  with  a  comma,  including  the  last  field. 

11.  Type  a  (s^  .  The  computer  should  now  be  running. 

If  a  field  is  simply  entered  as  a  comma  with  no  prefixed  number,  or  function  name, 
the  value  from  the  preceeding  test  will  be  used.  In  the  case  of  the  first  test,  do  not 
trust  this  value  to  be  zero. 

Function  Arguments 

NOT: 

1.  Selecting  this  function  causes  the  antenna  to  remain  stationary  at  the  given 
bias  angle.  No  input  parameter  is  used,  but  a  single  comma  must  be  typed  in  the 
parameter  field  (see  Fig.  1). 

POLY: 

1.  The  nth  field  in  the  argument  list  for  a  polynomial  function  is  the  coef¬ 
ficient  a  I  in: 

n- 1 

2  3  4  5 

a.t  a^t  aot  a.t  art 

^  ^  To'  100  1000  10000  100000 

where  t  is  given  as  elapsed  time  in  seconds  and  the  A  angle  is  in  degrees. 
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SINE: 

1.  Frequency  of  oscillation  in  cycles/sec.  (f). 

2.  Amplitude  of  oscillation  in  degrees  (A). 

A  )^  =  A  sin  (27rft) 

STEP: 

1.  Time  for  first  step  (seconds). 

2.  Height  of  steps  (degrees). 

3.  Repetition  time  (seconds). 

IMPU  LSE: 

1.  Time  for  first  impulse  leading  edge  (seconds). 

2.  Height  of  impulses  (degrees). 

3.  Duration  of  impulse  in  seconds. 

4.  Repetition  time,  leading  edge  to  next  leading  edge  (seconds). 

B.  Program  II  -  IBM  7094  Print  Program 

Tliis  program  reads  the  files  of  data  from  the  490  output  tape  and  prepares  a  BCD 
tape  (for  printing  on  the  1401).  For  each  file  on  the  490  tape,  a  print  file  is  generated 
on  the  BCD  tape.  Each  file  contains  an  identifying  label,  in  addition  to  the  edited  data 
obtained  from  the  490  and  converted  to  decimal  numbers  of  degrees. 

F'igure  4  shows  a  listing  of  a  control  deck  that  was  used  for  printing  the  output 
tape  produced  by  the  console  communications  of  Fig.  1.  The  first  card  gives  the 
number  of  files  to  be  listed  in  columns  9  and  10.  Tlie  seven  remaining  cards  (one  for 
each  file)  give  the  editing  factor  for  the  corresponding  files.  If  this  number  is  n,  every 
nth  point  on  tape  will  be  printed. 

Figure  5  shows  a  page  of  the  output  produced  by  the  specification  of  Fig.  4  and  the 
tape  from  the  490. 

Tlie  first  line  is  simply  the  run  label  that  was  typed  in  on  the  490.  The  columns 
are  (from  the  left): 

1.  Time  elapsed  in  seconds  from  first  application  of  drive  function  to  servos. 

2.  Azimuth  (in  degrees)  output  to  servos. 

3.  Elevation  output  to  servos. 

4.  Azimuth  input  from  encoders. 

5.  Elevation  input  from  encoders. 

6.  Difference  of  columns  2  and  4. 

7.  Difference  of  columns  3  and  5. 
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For  brevity,  only  a  single  page  of  output  is  reproduced  here;  however,  this  would 
normally  be  quite  copious. 

This  IBM  7094  program  is  a  relocatable  column  binary  deck  that  must  be  run 
under  FMS  control.  The  deck  is  labeled  "Haystack  Test  Print.  "  The  490  output  is 
mounted  on  A7.  The  print  tape  is  A6.  The  control  data  deck  is  placed  behind  the 
binary  deck  for  the  run. 

Control  Deck: 

1st  Card:  Col.  1  =  *,  Cols.  7-10  =  DATA 

2nd  Card:  Number  of  files  to  process  in  Cols.  1-10  as  a  right  justified 

integer. 

Rest  of  Deck:  One  card  per  file  to  process,  giving,  in  Cols.  1-10,  the 
frequency  of  selecting  points  for  printing  as  a  right  justified  decimal  integer  (1  = 
every  point,  2  =  every  2nd  point,  etc. ). 

C.  Program  III  -  IBM  7094  Plot  Program 

Tliis  program  reads  files  of  data  from  the  490  tape  as  they  are  addressed  by 
their  file  labels.  It  prepares,  as  output,  a  BCD  tape  for  driving  the  Calcomp  plotter 
(via  the  1401).  The  format  of  the  graphs  produced  is  quite  flexible. 

Figure  6  shows  a  listing  of  a  control  deck  that  was  used  to  plot  some  of  the 
information  from  the  output  tape  produced  by  the  console  communications  of  Fig.  1. 
The  first  card  gives  the  number  of  plots  to  be  made  in  Cols.  9  and  10.  Each  plot  is 
specified  by  a  pair  of  cards.  The  first  card  of  each  pair  contains,  in  Cols.  1-30, 
the  label  of  the  file  from  which  the  plot  is  to  be  made.  The  second  card  of  each  pair 
contains  six  variables,  in  fields  ending  in  Cols.  10,  20,  30,  40,  50  and  60,  respec¬ 
tively.  Tlie  first  of  these  provides  selection  of  the  coordinate  to  be  plotted  (0  means 
plot  azimuth;  1,  plot  elevation).  The  second  field  specifies  the  time,  in  seconds, 
between  plotted  points,  while  the  third  and  fourth  fields  state,  respectively,  the 
lower  and  upper  limits  of  the  time  axis  in  seconds.  The  fifth  field  gives  the  physical 
length  of  the  graph  in  inches  (measured  along  time  axis).  The  sixth  field  specifies 
what  data  is  to  be  plotted  for  the  given  file  and  coordinate.  This  appears  as  a  one¬ 
digit  code.  For  more  detailed  discussion,  see  p.  10. 

Figures  7  through  14  are  some  graphs  produced  by  the  data  cards  in  Fig.  6.  Note 
that  the  file  label  is  neatly  plotted  above  the  graph  and  that  the  fixed  length  of  the 
vertical  (angle)  axis  (5  inches)  makes  the  graph  suitable  for  full-size  reproduction  on 
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standard  8-1/2’'  x  11"  paper.  All  parameters  of  the  graphs  can  be  seen  to  conform 
to  the  relevant  specifications  of  Fig.  6.  Since  the  antenna  was  not  connected  to  the 
computer  at  the  time  these  tests  were  made,  Figs.  8,  JO,  and  14  show  errors  due  to 
delays  in  the  in-out  system,  but  do  not  show  true  servo  errors. 

This  IBM  7094  program  exists  as  a  relocatable  column  binary  declc  that  must  be 
run  under  FMS  control.  The  declc  is  labeled  "Haystack  Test  Plot.  "  The  490  output 
is  mounted  on  A7.  The  plot  tape  on  A6.  The  control  deck  is  placed  behind  the  binary 
deck  for  the  run. 

Control  Deck: 

1st  Card:  Col.  1  =  *,  Cols.  7-10  =  DATA 

2nd  Card:  Cols.  1-10,  right  justified  decimal  integer,  number  of  graphs  to 
be  plotted,  same  as  number  of  pairs  of  cards  which  follow. 

Rest  of  deck  consists  of  a  pair  of  cards  for  each  graph  to  be  drawn.  Of  each 

pair: 

2n  +  1st  Card:  Cols.  1-30,  left  justified,  the  exact  label  given  to  the  file  to 
be  plotted  when  that  file  was  generated  on  the  Univac  490. 

2n  +  2nd  Card:  Col.  10  =  0,  plot  azimuth;  1,  plot  elevation. 

Cols.  11-20  =  time  interval  between  data  points  to  be  plotted,  in  seconds,  as 
any  decimal  number,  right  justified. 

Cols.  21-30  =  lower  time  limit  of  plot,  in  seconds,  as  any  decimal  number, 
right  justified. 

Cols.  31-40  =  upper  time  limit  of  plot,  in  seconds,  same  format. 

Cols.  41-50  =  length  of  time  axis  on  graph,  in  indies,  same  format. 

Cols.  60  =  specifies  what  combinations  of  data  to  plot,  with  following 
interpretation: 

0,  nothing  plotted 

1,  only  input 

2,  only  output 

3,  input  and  output 

4,  only  error 

5,  input  and  error 

6,  output  and  error 

7,  input,  output,  and  error 
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IV.  INTERNAL  DETAILS  OF  THE  UNIVAC  490  ANTENNA  SERVO  TESTING  PROGRAM 


The  three  programs  will  again  be  discussed  separately.  This  material  should  be 
read  with  reference  to  the  program  listings  in  the  Appendices. 

A.  Program  I  -  Univac  490  Testing  Program 

This  program  supplies  azimuth  and  elevation  coordinates  to  the  antenna  servo 
every  4  ms  in  accordance  with  input  specifications.  The  time-dependent  variation 
of  the  coordinates  from  a  constant  bias  is  calculated  through  a  subroutine  for  the 
function  used.  The  manner  of  linking  subroutines  to  the  program  will  be  described. 

Execution  of  the  program  begins  at  ANTENTRY.  An  initial  remark  is  typed. 

By  use  of  the  subroutine  SINEFUNCS  a  table  of  the  sine  function  containing  2048  values 
is  calculated  and  stored  at  location  SINTBL.  This  corresponds  to  a  density  of  512 
points  per  quadrant.  The  table  is  subsequently  used  by  AZSINE  and  ELSINE  for 
driving  the  antenna. 

Next  an  interrogative  remark  is  typed  out  requesting  the  parameters  of  on-line 
printer  output.  The  first  ciaracter  typed  in  is  examined  to  determine  if  it  is  a  "Y" 
or  an  "N".  If  it  is  an  N,  the  program  jumps  to  INTCLKTEST  after  clearing 
ZXNLINEIND.  If  it  is  a  Y,  the  5th,  6th  and  7th  characters  are  entered  as  a  decimal 
number,  converted  to  octal,  and  stored  in  the  lower  half  of  ZXNLINEIND. 

The  entry  at  PROGRAM  is  used  to  write  an  end-of-file  mark  on  tape  #1  on 
completion  of  an  operation.  INPUTBUF  is  cleared  to  receive  the  run  label  which  is 
accepted  from  the  console.  The  label  is  converted  to  BCD  from  fieldata  by  FDTOBCD 
and  is  written  on  tape  #1.  The  program  now  accepts  the  running  parameters. 

The  line  of  data  is  first  tested  for  a  double  period  code  .  "  which  indicates 
termination  and  causes  a  jump  to  NOMORE.  If  this  is  not  present,  the  line  is  scanned 
from  left  to  right  and  separated  at  the  commas.  The  fields  are  stored  successively 
in  the  areas  designated  by  the  low-order  halves  of  TABLE.  The  operation  is  termi¬ 
nated  by  sensing  a  stop  "  ^  and  control  passes  to  VIRGULE.  The  program  calls 
DECON  to  convert  the  first  5  variables  listed  in  TABLE  to  octal  with  binary  points 
indicated  by  the  upper  halves  of  the  corresponding  words  in  TABLE.  If  a  variable  is 
indicated  only  by  a  comma,  the  number  in  the  register  is  left  unchanged;  otherwise, 
the  converted  number  is  placed  in  the  appropriate  location  in  CONTBL.  When  this 
is  complete,  control  passes  to  FUNIDLOOP. 

FUNIDLOOP  compares  the  bit  patterns  stored  in  FAZI  and  FELE  (addressed 
through  TABLE)  with  the  various  possible  names  in  FUNCTBL.  If  the  function  is 
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successfully  identified,  its  address  and  argument  count  is  transferrred  from  the 
appropriate  location  in  FUNADDTBL  to  FUNLOC  or  FUNLOC+1.  If  the  function 
name  cannot  be  identified,  an  appropriate  error  remark  is  typed  out  and  control 
is  returned  to  CONINETC. 

If  no  difficulties  are  encountered,  the  function  arguments  are  converted  by 
scanning  the  tables  appended  to  the  beginning  of  each  subroutine.  The  upper  half 
of  each  word  contains  the  location  of  the  binary  point  and  the  lower  half  the  address 
for  storage  of  the  converted  variable.  If  the  total  field  count  at  tlie  end  of  this  oper¬ 
ation  does  not  match  the  argument  counts  obtained  from  the  upper  halves  of  FUNLOC 
and  FUNLOC+l,  an  appropriate  remark  is  typed  and  control  passes  to  CONINETC. 
Otherwise,  the  real  time  section  of  the  program  is  initiated  at  ZKFLDCNT. 

AZST  and  ELST  are  set  equal  to  the  biases  of  azimuth  and  elevation  in  degrees, 
and  these  biases  are  converted  to  revolutions  with  a  binary  point  at  19,  The  time 
periods  specified  in  the  input  line  are  converted  to  a  count  of  servo  I-O  cycles  (4  ms 
each).  Tliey  go  into  NWAIT,  NRUN,  and  FREQOUT,  respectively.  A  loop  on  B1  for 
3750  cycles  (15  seconds)  moves  the  servos  to  the  initial  position  specified  by  the  input 
string. 

Tlie  drive  loop  is  initialized  with  two  index  registers  (B1  and  B4)  being  used  in 
tandem  to  count  running  times  longer  than  11111  ^  cycles.  POINTGEN  begins  the 
drive  loop.  It  ends  roughly  at  PRINTRET4-5.  In  this  loop,  the  elapsed  time  in 
seconds  is  calculated  and  stored  in  THYME;  the  appropriate  azimuth  and  elevation 
functions  are  called  and  the  results  (in  degrees)  stored  in  CONTBL.  These  quantities 
are  converted  to  revolutions  and  stored  in  AZIMUTH  and  ELEVATION.  Finally,  the 
computed  point  is  transferred  to  the  servos  via  DRIVEREAL  and  index  registers  are 
tested  to  determine  appropriate  action.  Register  B5  delays  output  until  the  servos 
can  start  to  track  the  applied  function;  i.  e.  ,  until  the  servo  response  has  reached  a 
steady-state  condition.  Register  B2  causes  a  line  of  output  to  be  added  to  the  output 
buffer  (through  WRITEBUF),  with  the  period  determined  by  FREQOUT  loaded  into  B2. 
B3  performs  a  similar  function  for  the  on-line  printer  (transferring  control  to 
PRINTER).  As  mentioned  before,  B1  and  B4  are  used  together  to  control  the  total 
number  of  cycles  of  the  loop. 

When  the  loop  is  finished,  whatever  remains  in  a  tape  output  buffer  is  written 
out  along  with  a  terminal  end-of-file  mark.  The  biases  in  CONTBL  and  CONTBLTl 
are  reset  to  their  initial  value  before  the  run. 
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Control  is  transferred  to  WRITEBUF  if  an  entry  is  to  be  made  in  a  tape  output 
buffer.  The  output  values  are  obtained  from  the  encoders  via  ENCODEREAL  and  the 
five  variables  TIME,  AZIMUTH,  AZOUT,  ELEVATION,  and  ELOUT  are  assembled 
into  six  words  so  that  a  blank  block  of  six  bits  will  appear  between  each  word  written 
on  tape  (this  extra  byte  makes  direct  reading  of  the  tape  on  the  7094  possible).  The 
data  is  assembled  into  BUFB  or  BUFA  depending  on  the  status  of  the  indicator  FIRST. 
When  B6  indicates  that  100  lines  have  been  entered,  TAPEWRITE  is  called  to  initiate 
the  output  buffer. 

At  the  end  of  a  tape  when  a  double  period  has  been  typed  in  as  an  input  line, 
control  passes  to  NOMORE  which  asks  if  continuation  is  desired.  A  YES  or  a  NO 
(or  a  Y  or  N)  may  be  typed  in.  The  program  also  asks  if  the  output  tape  should  be 
rewound  or  left  positioned.  Again,  YES  or  NO  will  be  accepted.  Tlie  indicated 
operation  is  performed  (rewind  with  interlock)  and  control  is  returned  either  to  MAIN 
or  TOPS  (137). 

Subroutine  TAPEWRITE  sets  up  the  channel  13  interrupt  location  (external)  and 
activates  the  appropriate  buffer  to  write  output  in  binary  high  density  mode  on  servo 
#1. 

Subroutine  DEGCON  divides  the  contents  of  the  A  and  the  Q  separately  by  360 
to  convert  degrees  to  revolutions  at  a  binary  point  of  19. 

Subroutines  DRIVE  and  ENCODE  are  dummy  routines  that  function  through 
jumps  to  DRIVEREAL  and  ENCODEREAL. 

TABLE  is  a  list  of  the  addresses  of  tlie  areas  for  storing  the  various  input  fields 
separated  from  the  input  line.  FUNCTBL  is  a  list  of  the  names  of  allowed  drive 
functions,  while  FUNADDTBL  contains  the  addresses  of  the  functions  and  the  length 
of  the  argument  storage  list  associated  with  each  one. 

ELPOLY  and  AZPOLY  are  the  polynomial  generating  subroutines.  Tliey  are  each 
prefixed  witli  a  list  of  the  addresses  of  the  six  coefficients  of  the  polynomials: 


0  -  0  = 
o 


Z 

n=0 


(p  -  (p^ 


=  Z 
n=0 


6+n 


where  f  =  (THYME/ 10),  0  =  azimuth  in  degrees,  and  cp  =  elevation  in  degrees.  The 
variable  £  is  stored  in  THYMEX.  0^  and  cp^  are  the  initial  biases  stored  in  AZST 
and  ELST. 
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Subroutine  DECON  converts  a  12-character  decimal  numeric  field  stored  in 
NUMBER  to  binary  according  to  the  placement  of  the  binary  point  specified  in  BINPNT. 
The  result  is  left  in  ACCUM.  Acceptable  forms  are: 

3. 1415926,  3,  3.  , 

+3.  1415926,  +3,  +3.  , 

-3. 1415926,  -3,  -3. , 

No  more  than  nine  numbers  may  be  used  in  the  input  field  and  it  must  be  terminated 
by  a  comma.  BINPNT  may  be  greater  than  31  or  negative,  if  desired;  however,  it  is 
usually  positive  and  less  than  31.  This  is  a  general  utility  subroutine  and  is  trans¬ 
parent  to  all  external  and  active  registers. 

FIXTC  is  the  error  recovery  program  for  magnetic  tape  interrupts. 

FDTOBCD  converts  a  word  in  the  Q-register  from  fieldata  to  BCD  by  use  of  the 
table  BCDTBL.  One  character  at  a  time  is  entered  into  B2  and  the  appropriate  code 
obtained  from  the  table  and  appended  to  ANZ. 

SINEFUNCS  computes  the  sine  of  an  angle  specified  in  revolutions,  with  sign 
corrected  in  all  quadrants.  Tlie  sines  are  actually  computed  in  the  first  quadrant  by 
the  Hastings  approximation  subroutine  SINE  FUN.  SINEFUN  obtains  the  sine  of  the 
fractional  revolution  stored  in  the  Q-register  at  B21.  The  result  is  left  in  the  Q  at 
B21. 

Subroutines  AZSTEP  and  ELSTEP  compute  the  step  functions.  This  is  accomplished 
by  finding  the  remainder  when  THYME-AZ(EL)  TIME  is  divided  by  AZ(EL)PERIOD.  If 
this  is  less  than  777g,  CONTBL(+l)  is  incremented  by  AZ(EL)HEIGHT. 

Subroutines  AZIMPULSE  and  ELIMPULSE  operate  similarly;  however,  both 
THYME-AZ(EL)TIME  and  THYME-AZ(EL)TIME  -  AZ(EL)DURAT  are  computed  and 
divided  by  AZ(EL)PERIOD.  If  the  former  remainder  is  sufficiently  small,  CONTBL(+l) 
is  incremented  by  AZ(EL)HEIGHT.  If  the  latter  is  small,  CONTBL(+l)  is  decremented. 

Subroutines  AZSINE  and  ELSINE  perform  a  table  look-up  operation  to  simulate  the 
sine  without  actual  computation  with  the  Hastings  approximation.  The  product 
(AZ(EL)OMEGA*THYME)  is  calculated,  modulo  1  at  B21.  The  high  order  bits  in  B1 
are  used  to  address  SINTBL,  while  the  low  order  11  bits  are  used  to  interpolate  between 
successive  table  entries.  The  result  is  multiplied  by  AZ(EL)AMPLITUD  and  added 
to  AZ(EL)ST,  in  the  A-register. 

Subroutines  AZ RANDOM  and  ELRANDOM  provide  for  pointing  the  antenna  in  a 
fixed  direction.  They  are  referenced  whenever  the  driving  specifications  AZNOT  and 
ELNOT  have  been  entered. 
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Subroutine  PRINTER  effectively  converts  the  azimuth  and  elevation  commands, 
responses,  and  errors  to  decimal  or  octal  and  causes  them  to  be  printed  on-line. 

The  conversion  of  a  line  of  data  containing  7  quantities  (decimal)  or  9  quantities 
(octal  mode)  is  multiplexed  with  the  main  pointing  calculations.  Tlie  list  of  JPTABLE 
specifies  the  order  of  execution  of  the  various  operations  for  decimal  mode  printing. 
Each  time  PRINTER  is  entered,  control  passes  to  the  next  subroutine  in  the  list. 
GETNUMS  picks  up  the  current  values  of  the  parameters  to  be  printed  and  stores 
them  in  FIELDS.  FDVAR  is  used  to  transmit  information  from  the  decimal  and  octal 
conversion  subroutines.  The  fieldata  line  to  be  printed  is  assembled  in  PRBUFER. 

Each  section  named - CON  (example,  TIMECON)  converts  the  current  value  of 

the  integer  part  of  the  specified  variable  to  decimal  fieldata.  Each  section  named 

- CONB  (example,  TIMECONB)  performs  the  binary-to-decimal  conversion  to  four 

decimal  places  for  the  fractional  part  of  the  same  variable.  The  result  is,  in  each 
case,  stored  in  the  appropriate  area  of  PRBUFER.  RUNPRINTER  causes  the  contents 
of  PRBUFER  and  STOPMARK  to  be  transmitted  to  the  on-line  printer  control  system. 

It  also  terminates  the  multiplex  routine  PRINTER  by  clearing  the  pointer  in  PRINTINDIC. 

PRCONVl  makes  all  numbers  received  through  location  NUMBER  positive,  splits 
them  into  integer  and  fractional  parts,  and  calls  INTEGERCON,  which  converts  the 
integer  part  to  a  signed  three-digit  fieldata-decimal  integer  terminated  by  a  decimal 
point.  FRACTCON  converts  the  fractional  part  stored  in  FRACTION. 

DRIVINIT  sets  up  interrupt  locations  for  the  azimuth  and  elevation  output  buffers. 
DRIVEREAL  transmits  data  to  the  azimuth  and  elevation  interface  systems. 

Tlie  following  table  gives  the  storage  locations,  binary  scale  factors,  and  units 
for  various  quantities  used  frequently  in  the  testing  program. 

Quantity 

Azimuth  bias,  degrees,  B18 
Elevation  bias,  degrees,  B18 
Azimuth,  degrees,  B18 
Elevation,  degrees,  B18 

Azimuth,  revolutions,  B19  (command  number  system) 

Elevation,  revolutions,  B19  (command  number  system) 

Azimuth  encoder  bit  pattern  (command  number  system) 


Location 

AZST 

ELST 

CONTBL 

CONTBL+1 

AZIMUTH 

ELEVATION 


AZOUT 

Elevation  encoder  bit  pattern  (command  number  system)  ELOUT 
Relative  real  time,  seconds,  B18  THYME 
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ENCODEREAL  reads  the  azimuth  and  elevation  encoders  from  the  interface  equipment 
into  locations  AZOUT  and  ELOUT,  performing  a  small  calculation  to  correct  for  angles 
in  the  overlap  zones  in  azimuth.  Subroutine  DRIVFIX  answers  interrupts  associated 
with  the  azimuth  and  elevation  channels.  Subroutine  TYPEIN  answers  interrupts 
associated  with  the  console  keyboard  and  performs  certain  elementary  operations.  As 
an  example,  if  O  or  D  is  typed,  it  transfers  control  to  FIXOCTPRNT  or  FIXDECPRNT, 
respectively,  which  set  the  printing  mode  to  octal  or  decimal,  as  indicated.  It  performs 
a  similar  function  with  respect  to  the  characters  S  and  P  and  the  locations  STOPRINT 
and  GOPRINT. 

Subroutine  RUNNY  is  a  substitute  routine  for  RUNPRINTER  in  the  sense  that  it 
terminates  the  print  conversion  computation  without  printing  a  line.  JPTABLEOCT 
performs  the  same  function  in  octal  printing  as  does  JPTABLE  in  decimal  printing. 

Each  routine  named - OCT  calls  CONOCT  to  convert  the  prefixed  variable  to  fieldata 

octal  (ten  characters).  As  an  example,  AZDIFOCT  converts  the  error  in  azimuth  to 
octal  for  printing. 

Since  the  fields  for  octal  and  decimal  printing  do  not  precisely  coincide,  BUFKLEER 
ensures  that  extraneous  characters  are  not  carried  over  from  one  format  to  the  other 
if  format  is  changed  in  mid-run. 

Addition  of  Drive  Function  Subroutines  to  the  Testing  Program 

A  drive  function  subroutine  may  be  easily  added  to  the  program  by  preparing  it  in 
the  proper  format  and  by  making  a  few  entries  in  some  tables. 

Tlie  call  name  (for  console  type-in)  for  the  subroutine  should  be  entered  in  fieldata 
in  FUNCTBL.  Tlie  length  of  the  argument  list  and  the  entry  address  must  be  entered  in 
FUNADDTBL. 

Each  subroutine  must  be  prefaced  with  a  table  of  the  addresses  and  binary  points 
of  its  input  arguments  (excluding  the  time  variable).  The  binary  point  is  placed  in  the 
upper  half-word  and  the  variable  address  in  the  lower.  The  first  argument  typed  in 
will  be  placed  in  the  last  specified  address,  etc.  ,  so  that  the  variable  addresses  are 
listed  in  the  table  in  the  reverse  order  from  the  typed-in  line. 

B.  Program  II  -  IBM  7094  Print  Program 

This  program  converts  the  490  output  tape  to  a  BCD  tape  that  will  print  on  the 
1401.  The  number  of  files  to  be  processed  from  input  tape  A6  is  read  from  A2  and  a 
DO- LOOP  is  initialized  for  this  number  of  iterations. 
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For  each  file  to  be  printed,  the  file  label  is  read  from  A7  and  written  in  BCD 
mode  on  A6.  The  editing  factor,  N,  is  read  from  A2.  From  this  point  the  program 
reads  in  490  records  of  500  36 -bit  words,  converts  them  and  writes  them  out  until 
an  end-of-file  is  encountered  on  A6,  whereupon  the  output  file  on  A7  is  terminated  and 
the  index  of  the  major  DO  incremented  (statement  101). 

If  an  end-of-file  is  not  encountered  by  subroutine  READER  in  reading  A6,  lOF 
will  be  zero  and  the  program  will  reverse  the  order  of  the  input  array,  storing  the 
reversed  array  in  BUFFER.  The  entire  arr^y  is  converted  to  floating  point  by  sub¬ 
routine  FLOTER,  and  the  various  terms  in  the  array  are  scaled  by  the  appropriate 
factors  (.004  for  time,  .000686645507  for  angles).  The  errors  in  azimuth  and 
elevation  are  computed  and  stored  in  DIFFER.  Finally,  every  Nth  point  is  written 
out  on  A7  in  an  appropriate  format  for  1401  printing. 

C.  Program  111  -  IBM  7094  Plot  Program 

This  program  converts  the  490  output  tape  to  a  BCD  tape  suitable  for  driving  the 
Calcomp  plotter  to  produce  graphical  output  in  accordance  with  input  specifications. 
First,  the  program  requests  mounting  of  an  input  tape  on  A7  and  pauses.  When 
restored,  it  reads  the  number  of  graphs  to  be  plotted  from  A2,  initializes  the 
Calcomp  subroutine,  PLOTS,  and  begins  a  DO-LOOP  for  the  number  of  frames  of 
three  graphs  each  to  be  plotted. 

For  each  graph  to  be  plotted,  A7  is  rewound.  Tlie  relevant  file  label  is  read 
from  A2,  along  with  the  editing  parameters.  Subroutine  READER  is  called  to  read 
the  490  output  tape  (A7),  placing  a  file  label  in  BUFXIN.  The  order  of  the  words  is 
reversed  as  they  are  transferred  to  BUFRIN.  The  file  label  from  A7  is  compared 
with  the  desired  label  from  A7.  If  they  do  not  match,  A7  is  advanced  one  file  by 
FILS  PA  and  the  process  repeated. 

When  the  relevant  file  has  been  located  in  this  manner,  control  passes  to  state¬ 
ment  24.  The  program  section  from  24  to  211  is  identical  to  the  conversion  section 
of  the  program  CONVERT  described  above.  After  a  record  has  been  converted,  the 
time  fields  are  scanned  to  find  one  greater  than  TIML.  When  this  is  found,  the 
appropriate  line  of  data  is  transferred  to  a  plot  buffer,  and  TIML  is  advanced  by 
SPACE.  The  process  continues  until  the  input  record  is  exhausted,  then  another  is 
read  (statement  240).  Eventually,  TIML  exceeds  TIMU,  and  the  remainder  of  the 
file  is  spaced  over. 
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Since  the  Calcomp  scaling  routine  requires  an  array  of  non-identical  entries, 
the  plot  buffers  are  checked  for  this  before  calling  SCALE.  Tliose  arrays  which  are 
scalable  and  plotable  are  processed  in  accordance  with  indicator  lUJ.  After  the  lines 
of  data  have  been  plotted,  AXIS  and  XAXIS  are  called  to  draw  in  the  relevant  axes. 
Finally,  the  label  is  plotted  by  SYMBL4,  and  the  graph  origin  is  reset  for  the  next 
plot. 

When  three  graphs  have  been  drawn,  the  paper  is  advanced  three  inches  beyond 
the  end  of  the  longest  graph  and  the  pen  is  reset  to  the  right-hand  margin.  After  all 
the  graphs  specified  by  the  input  deck  have  been  plotted,  the  output  tapes  are  rewound 
and  a  notice  to  dismount  is  printed  on-line. 

Prior  to  the  final  exit,  the  program  pauses  (HPR77777). 
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APPENDIX  B 
7094  PRINT  PROGRAM 


*  LISTS 
CCONVERT 

C  this  program  reads  a  tape  prepared  on  the  UNIVAC  *90  AMD  PRINTS  IT 

C  WITH  SOME  SMALL  CALCULATIONS 

DIMENSION  BUFFERISOO  ) . DI FFER ( 200 ) . BUFFRX ( 500 ) 

READ  INPOT  TAPE  2»100fNTlM  _ 

100  FORMAT  (IlOI  _  _ 

00  101  III-ltHTIM  -  -  ^ 

10'  t0F«0 - 

CALL  READER  ( BUFFRX ( 5 ) *5  * lOF ) 

DO  102  NN«1,5 

- - ^ — 

lOZ  BUFFER  (RN}«BUFFRX(MHHN}^ - 

WRITE  OUTPUT  TAPE  6*20*  (  BUFFER  jL  J  »).«i  t5  ) 

EO  FORMAT  (5A6) 

READ  INPUT  TAPE  2»6»N 
6  FORMAT  (110) 

_ H«o  -  ^ - 

t?  CALL  READER  (BUFFRXI500I  iSOOtlOF) _ ^ - 

IF  IlOF)  11»10*»11  -  = - 

10*  DO  103  NNi=l*500  - 

NNHN»501-NN 

103  BUFFER(NN)«BUFFRX(NNNN) 

-  De-2"LOEX*l»50O  —  = — 

E  CALE  FLOTER  ( BUFFER  lEOEX  H  =  — - 

oo  3  LOEX-1  ♦  *96 . 5 - 

BUFF^  ( LOEX ) -buffer  ( LOEX  t  *  .  00  * 

DO  *  NDEX=1»* 

JDEX-LDEX+NDEX 

r  *  BDFFERIJDEXl-BOFFERf  J0EX)»iSeD«66*S507 _ 

"  ^  gOBTIMUE  —  _  - 

-=-m  5  EDEX»T»100  _ ! _ 

JDEX«2*LOEX-t 

NDEX«5*LDEX-3 

MDEX-NDEX+1 

DIFFER ( JDEX ) -BUFFER { NDEX ) -BUFFER ( MDEX ) 

5  DIFFER(JOEX*ll-BUFFERtNOeX4-2I-BUFFER(MDEX*2)  _ 

— OQr?  LOEX-ltlOO -  : -  “ 

— - -TlOEX«XMODF(M«Nr -  ^ - 

M-M  +  1 

IF  (JDEX)  7»8»7 

8  NOEX-5»LOEX-* _ _ 3_ 

=^=^iDEx*2»fceEX!*i  -  -  ■  — ' — 

I XOEX- lOEX^l  ■  “ 

WRITE  OUTPUT  TAPE  6 *9 ♦ ( ( BUFFER ( T ) . I -NDEX .NXDEX ) * ( DIFFER ( I ) ♦ I « IDEX * 
IIXDEX) ) 

9  FORMAT  (9X#5F15.9*2Et5»8»  - 

^ - T^COWTIMUE  -  =  - 

11  END  FTLE=T —  ^  ~  ~ 

101  CONTINUE 
CALL  EXIT 

END  ^ _ 

*  FAP  — ^ 
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APPENDIX  C 
7094  PLOT  PROGRAM 


LISTS 
=CPlOTTFR 

C 
C 


THIS  PROGRAM  ACCEPTS  THE  CONVERTED  UNIVAC  490  TAPE  AND  PLOTS  THE 
DATA  IN  ACCORDANCE  WITH  THE  INPUT  SPECIFICATIONS. 

OTMENSION  BUFFER  t  JOOO 1 1 BUFRIN ( 500 1 . BUF* I N ( 500 1 . DIFFER ( 200 1  ♦” 

RL0T81 1 5000 J  »PL0TB2( 50 


Z0©^^^B3(5OOO)  .PLOTTOFlgOirT»XLARFL  1 5 1 .XLARSR (  5 )  .XLFN  (  5 1 


PftINT=I^=  - -  = - 

WRITE  OUTPUT  TAPE  3*15 

15  FORMAT  {46H1M0UNT  INPUT  TAPE  ON  A7  AND  SCRATCH  TAPE  ON  A6 ) 
PAUSE  70707  — — 

READ  INPUT  TAPE  2il»NGRAPH 
1  FORMAT  I 1101 
XSTA«0.0 
YLEN*6.0 

CALL  PLOTS  (BUFFER (1000) *1000) 


NGR-t<N0RAPH-l|/3I<-l 

DO  2  1«1.N6R 

CALL  PLOT  fO.0.“29.O.-3) 

CALL  PLOT  (0.0, 2. 5.-3) 

DO  3  J«l,3 
NTEST=( ( I-l )*3)+J 
rSTA«0.0 
REWIND  7 

READ  INPUT  TAPE  2 *21 , (XLABSR( L ) »L«1 ,5 ) 

READ  INPUT  TAPE  2 *4. JTESTiSPACE*TIML*TIMU*XLEN( J J ♦ lUJ 
Z1  FORMAT  (5A6) 

4  FORMAT  (I10*4F10.0»I10) 

25  CALL  READER  ( BUFXIN ( 5 ) *5 » lOF )  ■ 

DO  lOO  NN-1.5 

100  XLABELfWO-BUFXINlNNNN) 

DO  22  L-l*5 


S  XLABEL(L)»(XLABEL(L)*( -XLABSR(L) ) )+( XLABSR ( L ) *( -XLABEL ( L ) ) ) 

IF  tXLABEL(L»)  23*22,23 

22  CONTINUE 
GO  TO  24 

23  CALL  FIL5PA 
GO  TO  25 

24  KN«1 
240  INDIC*! 

CALL  reader  (BUFXIN (500) *300* IOF) 

IF  (IOF)  205*206*205 

206  DO  207  NN=1*500 
NNNN»501“NN 

207  BUFRIN(NN)«BUFXIN(MNNN ) 

DO  208  tDEX«l*500 

208  CALL  FLOTER  (BUFR1N(LDEX ) ) 

DO  209  L0EX*I*496*5 
BUFRIN(LDEX)»BUFR1N(LOEX)#.004 
DO  210  NDEX-1.4 


JOEX«NOEX+LDEX 

210  eUFRINl JDEX )=eUFRlN ( JDFX )*, 000686645507 

209  CONTINUE  ,  _ 

DO  21t  LDEX«1*100 
JDEX«2*LDEX-1 
N0EX«5*LDEX-3 
MDEX=NDEX+1 


i 
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OrFFER ( JDEX  J ■BUFRI N ( NDEX ) -BUFR I N ( MDEX ) 

2J1  DIFFER (JOEX+H«eUFft IN (Noex+2) -BUFRI N(MDEX+2)  ^ 

218  DO  212  LDEX«INDIC*496*5 
IF{BUFRIN(LDEX)-TIML)  212*213 *213 

212  CONTINUE 
60  TO  240 
=213  INDTC-LOEX 

■TlMl.«TlMt.+sPACE 
PlOTBl (KN)»BUFrin(INDIC) 

LDEX=<2*( INDIC-l)/5)+l 
IF(JTEST)  214.215*214 
21^  FL0TB2(KN)=BUFRIN( INDIC+3) 

PL0TB3(XN)»BUFRIN( INDrc+4l 
PL0TB4<KN)-0IFFER<LDEX+1 ) 

GO  TO  2l6 

215  PL0TB2(XN)*BUFRIN(INDIC+1) 

PL0TB3(ICN)=BUFRIN(  INDIC+2) 

PL0TB4 ( KN ) =DrFFER ( LDEX ) 

216  KN=KN+1 
INDIC-INOIC+5 

IF  (5000-XN)  217.219.219 

219  IF  (TIML-TIMU)  218,218*217 

217  CALL  FILSPA 
205  162*1 

KN»KN-1 

DO  100*  KXY=1*<N 

IF(PLOTB2(KXY)-PLOTB2(l) )  1005*1004,1005 

1004  CONTINUE 
IB2=0 

1005  IB4*1 

DO  10'06  KXY*1*KN 

IF  (PL0TB*(KXy)-PL0TB4(l))  1007*1006*1007 

1006  continue 

IB4«0 

1007  CALL  SCALE(PL0TB1*KN*XLEN(J)*XMIN*DX) 

IF<I82+IB4)  1014*12*1014 

1014  IF<IB2»  1010*1011*1010 

1010  CALL  SCALE<PL0TB2*KN*5.0*YMIN2*DY2 ) 

DO  16  K*1,XN 

16  PLOTB3(X)=(PLOTB3(K)-YMIN2)*10./(DY2) 

1011  IF(IB4)  1012*1013*1012 

1012  CALL  SCALE(PL0TB4,XN,5,0*YMIN4,DY4) 

1013  IF(rB2)  1040,1041*1040 

1040  IF  nuj)1101*12,llOl 

1101  IF  (XM0DF(I0J,2) )  1102*1103*1102 

1102  CALL  LINE(PL0TB1,PL0TB2,KN) 

1103  IF  (XMODF( IUJ,4)  -  2)  1041*1104*1104 

1104  CALL  LINE(PLOTBl*PLOTB3*KN) 

1041  IFIIB4)  1042*1043*1042 

1042  IF  (XMODF( IUJ*8)  -  4)  1043*1105*1105 

1105  CALL  LINE(PLOTB1*PLOTB4*KN) 

1043  CALL  AXIS(XSTA*YSTA*4HTIME*4*XLEN(J) ,OiO*XMIN*DX 
XSTAX=XSTA-0*5 

IF  (JTEST)  11*10*11 
10  IF(IB2)  1020*1021*1020 

1020  CALL  AXIS(XSTA»YSTA*7HAZIMUTH*7*5.0*90.0*YMtN2,DY2  ) 

1021  IF(IB4)  1022*12*1022 

1022  CALL  XAXIS(XSTAX*YSTA,13HAZIMUTH  ERROR  *  1 3  *  5 . 0  *  90 • 0  * YMI N4 *DY4 *  3 ) 


60  TO  12 
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iYSTA  f  9mSSVAti  Oti*9*5, 0*90.0*  Ym»ti6y  2  } 

==ietr^ffe4t^««m^to«^ -■  -  -- 

1032  CALt  XAXIS(XSTAX»YSTA*15HELEVATION  ERROR . 15 ♦ 5 .0» 90 .0 ♦ YM IN4 ♦DYA *3  ) 
12  X-XSTA+3.0 


r - tatg-matCr^iri 


CALL  SYMBLAtXtY*. 21, XLABEL(5) 40*0*30) 
CALL  PLOT  (0.0»8.5*-3) 


-  -  -tf  «l«Tg5T-H6HAPm  _ -  -- 

^ - 3 

■.- 

-  ^  WTA  J»li3 

IF  (IZ)  32,31,32 
32  IZ-0 

_ : _ = 

'  SB  TO  14  -  = 

«[tHitJt-)a.Efii4svn  t*»i**iT 
rrusv  -j 

lA  CONTINUE 

XLEN(JSV)-XLEN( JSV)+A,0 


-  ^  -CALL  T»tOmLCIITVSVl aO* 0t*3» 


13  PRINT  30 

WRITE  OUTPUT  TAPE  3,30 


-  aOUMBMAT  < A9H101 SWOUHT  TAPE  A6  AHB  PLOT  IT*  SAVE  INPUT  TAPE  AT) 
77777  
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APPENDIX  D 


_ ^DIFIED  AXIS  PLOTTING  PROGRAM 

- 

CXAXIS 

SUBROUTINE  XAXIS  (X ♦YtBCDtNC# S IZE » THETA »YMIN tDY t IPEN ) 

TH  «  THETA  /  57#29578 

N*SI2E‘»*0.50 

YB*SINF(TH1 

XA«x-o*n»yB 

XB*C05ErtH)  ‘ _ 

XC«X  - 

YA«Y-0*1*XB 

YC*Y 

CHAR«ABSF{Y«IN)  * 

VALUE* ABSF  <YHIN+OY ) 

IF(CHAR-VALUC)5*6#6 

5  CHAR-VALUE 

6  N1*0 

VALUEhIOOOO* 

13*3 

14  IF(CMAR-VALUE)15#16*16 

15  Nl-Nl+1 
VALUE*VALUE/10.0 
GO  TO  14 

16  DO  20  I»1*N 

CALL  PL0TfXA#YA*I3) 

CALL  PLOT(XC»YC»IPEN) 

XC*XC4XB 

YCkYC+YB 

XA-XA+XB 

YA»YA+yB 

CALL  PLOTfXC*YC#IPEN) 

20  I 3* I PEN 

- CALL  PLOTrXA,YAfIPEN) 

XA*XC~*12 

YA»YC-.12 

XC*XC*XB 

YC*YC*YB 

N*N+1 

M  30  I«1*N 

VALUE*  t tX C+ YC ) *0Y/ 1 0 • 0 ) + YM  t  N 

CALL  NUMBER (XA*YA fOt 104V ALUE* theta ♦Nl ) 

XA»XA~XB 
YA*YA-YB  _ 

XC*XC“XB 
50  TC*YC-YB 
VALUE=NC/2 

XC*X+XB*{SlZE/2t0-0#12*VALUE)-0t3*YB 
YC*Y+YB*(SlZE/2*0-0.12*VALUE)«0.46*XB 
CALL  SYMBL4<XC«YC*0*l4fBCD*THETAfNC) 

RETURN 

END 
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APPENDIX  E 

TAPE  READING  SUBROUTINE 


ENTRY 

READER 

READER 

CLA 

1»4 

$TA 

CWORD 

CLA* 

2i4 

STD 

CWORD 

RTBA 

7 

RCHA 

CWORD 

*  .’  • 

TCOA 

« 

TEFA 

END 

STZ* 

3»4 

TRA 

4i4 

END 

CLA 

•1 

STO* 

■3#4 

TRA 

4»4 

'  ■  CWORO 

tORP 

»  • 

lOCD 

END 

f  « 
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#  WAP 
-- j  ^  EHipr 
FtCTER  CAI* 

ALS 

PBT 

THA 


APPENDIX  F 

mxitiaxm  point  cxinversion  subroutine 


F^&T» 

*■^7 


_ 

Ctit 

*1 

Sf4 

2te 

CAL* 

TiA 

ORA 

•07' 

COM 

-  TilA 

*+2 

ri£z  ^ 

€AL* 

It* 

— — — 

■OZI 

-  . 

~5Lvr* 

^lA 

CLA 

ZI6 

TIE 

PXA 

FS8* 

ZA6 

*f 

1*4 

ZTC 

-  PXA 

710^ 

f  * 

FAO» 

1*4 

ZOG 

STO* 

1*A 

ZfS 

TRA 

PZW 

2*4 
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APPENDIX  G 


FILE  SPAaNG  SUBRCXJTINE- 


RCHA  CWORD 
TCOA  # 

TCFA 

TRA  *-4 

_ TRA  lf4 

CWORO  IWm  M 

ioco 

EHD 
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READY 

LDaS2a32 

READY 

PSai00/5^  ti 


-62-2863 


THIS  13  THE  AnjTENNA  TESTING  PR03RAN. 

DO  YOU  WISH  TO  HAVE  OUTPUT  ON  THE  ON-LINE  PRINTER. 
IF  SO,  HOW  OFTEN. 

YESn2  50  “o 

TYPE  TEST  IDENTIPIOATION. 

SINE  AND  STEP  pUNOTIOMS® 

TYPE  TEST  PARAMETERS. 


90  ,0  , 1  , 100  ,  1  , AZ3I VE,ELSTEP, .02,5,5,3,7.5,® 

TYPE  TEST  identification. 

POLY  AND  IMPULSE  p'J  NOT  IONS® 

TYPE  TEST  PARAMETERS. 

,4  5,  ,3  5,  ,  AZPOLY,ELIMP,0  ,  72  .6  ,-  46 . 49  52  1  t,  12  .9  ,  -  1  .65,  .  AS  , 
5,10, .5,10,® 

TYPE  TEST  IDENTIFICATION. 

CONSTANT  3IASES® 

TYPE  TEST  PARAMETERS. 

92 .1393  75 ,42.390525  ,,  10  ,1  00  ,AZN3T  ,  ELECT  ,  ,  ,  ® 

TYPE  TEST  IDENTIFICATION. 


END  OP  DEMONSTRATION® 

TYPE  TEST  PARAMETERS. 

. .  ® 

DO  YOU  WISH  TO  CONTINUE 
NO® 

DO  YOU  WISH  TO  REWIND  OUTPUT  TAPE. 

YES® 

READY 


Fig.  1.  Sample  console  communications. 
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-62-28641 


i  •  W9 

90*6241 

91*2407 

90*6241 

91*2380 

0*0000 

-0*0027 

0*0000 

0*0000 

0*0000 

0*0000 

O.OOOO 

0  *  OOOO 

0-0000 

0*0000 

0  *  0000 
0-0000 

2*9959 

3.9959 

9l *8381 
92*4060 

91  *8354 
92.4039 

— 0*  0027 
-0.0020 

0*0000 

0*0000 

4.9959 

92.9367 

92.9347 

-0.0020 

0.0000 

0.0000 

0-0000 

5*9959 

93*4208 

93*4188 

-0«0020 

2*9999 

2*9999 

0*0000 

9*9959 

7*9959 

93*8507 

94*2201 

93*8493 

94*2187 

-0*0013 

-0*0013 

2-9999 

2*9999 

2.9999 

2*9999 

2*9999 

2.9999 

0-0000 

0*0000 

0*0000 

0.0000 

9*9959 

9,9959 

10*9959 

11.9959 

12*9959 

13. *959 

94*5229 

94,7543 

94.9108 

94.9898 

94*9898 

94.9115 

94*3215 

94,7536 

94*9102 

94,9891 

94*9898 

94.9122 

-0*0013 

-0.0006 

-0.0006 

-0.0006 

o*'oooo 

0.0006 

2*9999 

2.9999 

2.9999 

2.9999 
5*9999 

5.9999 

2.9999 

2.9999 

5*9999 

5*9999 

0-0000 

0.0000 

0.0000 

0*0000 

14*9959 

94*7537 

94.7563 

0*0006 

5.9999 

5.9999 

0*0000 

15*9959 

94.3249 

94,5236 

0.0006 

5*9999 

5.9999 

0.0000 

16.9959 

17*9959 

94.2228 

93*8541 

94.2242 

93*8555 

0*0013 

O.OOI3 

5.9999 

5.9999 

5*9999 

5.9999 

0*0000 

0.0000 

18*9959 

93*4243 

93*4263 

0*0020 

3.9999 

5*9999 

0*0000 

19*9959 

92*9409 

92*9429 

0*0020 

5*9999 

5.9999 

8*9998 

8.9998 

0-0000 

20*9959 

21.9959 

92*4108 

91*8429 

92*4128 

91*8450 

0*0020 

0.0020 

8*9998 

8.9998 

0*0000 

0*0000 

22.9959 

91.2455 

91.2483 

0.0027 

8.9998 

8.9998 

0-0000 

23.9959 

90*6289 

90*6310 

0*0020 

8.9998 

8.9998 

0*0000 

24.9959 

90*0020 

90*0048 

0.0027 

8*9998 

8*9998 

0*0000 

25*9959 

89*3600 

89*3627 

0*0027 

8*9998 

8.9998 

0*0000 

26.9959 

88.7441 

88.7461 

0*0020 

8.9998 

8*9998 

0*0000^ 

27.9959 

88. 1474 

88* 1494 

0*0020 

11.9998 

11.9998 

0*0000 

28.9959 

87,5795 

87,5816 

0.0020 

11.9998 

11.9998 

0.0000 

29.9959 

87.0501 

87.0522 

0*0020 

11.9998 

11.9998 

0*0000 

30.9959 

86,5674 

86,5695 

0.0020 

11*9998 

11*9998 

0*0000 

31*9959 

86*1389 

86*1403 

0.0013 

11*9998 

11.9998 

0*0000 

32,9959 

85,7709 

85,7723 

0.0013 

11.9998 

11,9998 

0.0000 

33.9959 

85.4701 

85.4708 

0*0006 

11.9998 

11.9998 

0*0000 

34.9959 

85.2401 

85,2415 

0.0013 

11.9998 

11.9993 

0.0000 

35*9959 

83*0856 

85*0863 

0*0006 

14.9997 

14.9997 

0*0000 

36*9959 

85*0087 

85*0087 

0,0000 

14*9997 

14*9997 

0.0000 

37*9959 

85*0101 

85*0101 

0*0000 

14*9997 

14*9997 

0*0000 

38.9959 

85.0904 

85.0904 

0.0000 

14.9997 

14*9997 

0*0000 

0*0000 

39.9959 

85.2484 

85*2477 

-0.0006 

14.9997 

14*9997 

40*9959 

85*4811 

85.4798 

-0.0013 

14*9997 

14.9997 

0-0000 

41*9959 

85*7846 

85*7833 

-0.0013 

14.9997 

14*9997 

0-0000 

42*9959 

86*1554 

86. 1534 

-0*0020 

17*9997 

17*9997 

0 . 0000 

43*9959 

86*5859 

86*5846 

-0.0013 

17.9997 

17.9997 

0  *  0000 

44*9959 

87*0707 

87*0687 

-0*0020 

17*9997 

17*999? 

0*0000 

45*9959 

46*9959 

47,9959 

48*9959 

87.6022 

88.1707 

88,7688 

89.3834 

90,0123 

87.6001 

88.1687 

88,7660 

89.3833 

90.0096 

-0.0020 

-0.0020 

-0.0027 

-0.0020 

-0.0027 

17*9997 

17.9997 

17.9997 
17*9997 
17,9997 

17.9997 
17.999? 

17.9997 
17*9997 
17,9997 

0*0000 

0*0000 

0.0000 

0*0000 

0*0000 

49.9959 

50*9959 

90*6234 

90.6214 

-0*0020 

20*9996 

20*9996 

0*0000 

51.9959 

91,2407 

91,2380 

-0.0027 

20.9996 

20*9996 

0.0000 

52*9959 

53.9959 

91.8381 

92,4060 

91.8354 

92,4039 

-0.0027 

-0.0020 

20*9996 

20.9996 

20.9996 

20.9996 

0-0000 

0.0000 

54*9959 

92*9367 

92.934? 

-0.0020 

20*9996 

20*9996 

0*0000 

55*9959 

56*9959 

57*9959 

93.4208 

93.8507 

94.2201 

93,4188 

93*8493 

94.2187 

-0,0020 

-0*0013 

-0.0013 

20*9996 

20*9996 

23.9996 

20.9996 
20*9996 

23.9996 

0.0000 

0-0000 

0.0000 

58.9959 

94*5229 

94.5215 

-0.0013 

23.9996 

23*9996 

0-0000 

59*9959 
^  60*9959 

~  61*9959 
42*9959 
63*9959 

94*7543 

94*9108 

94*9898 

94*7336 

94*9102 

94*9891 

-0.0006 

-0*0006 

-0*0006 

23*9996 

23.9996 

23*9996 

23*9996 

23*9996 

23*9996 

0*0000 

0*0000 

o.oooo 

94*9898 

94.9115 

94.9898 

94*9122 

0*0000 

0.0006 

23*9996 

23*9996 

23*9996 

23.9996 

V  •  VWV 
0-0000 
0-0000 

64*9959 

94*7557 

94.7563 

0.0006 

23*9996 

23.9996 

0.0000 

65*9959 

66*9959 

-  94*5249 

-  94*2228 

94.5256 

94*2242 

0*0006 

0.0013 

26*9995 

26.9995 

26.9995 

26*9995 

0*0000 

0*0000 

67*9939 
_  ^.9959 
69*9959 

93*8541 

93*4243 

92*9409 

93*8555 

93*4263 

92.9429 

0-0013 

0*0020 

0.0020 

26-9995 

26*9995 

26.9995 

26.9993 

26.9993 

26.9995 

0*0000 

0-0000 

0*0000 

70*9959 

92.4108 

92.4128 

0*0020 

26,9995 

26,9995 

0.0000 

71*9959 

^72.9959 

91.8429 

91.2435 

91,8450 

^91,2483 

0*0020 

0*0027 

26*9995 

29,9995 

26.9993 

29.9995 

0*0000 

0.0000 

73*9959 

74.9959 

75*9959 

90*6289 

90.0020 

89,3600 

90.6310 

90.0048 

89.3627 

0-0020 

0*0027 

0.0027 

29*9995 

29*9995 

29.9995 

29*9995 

29.9995 

29.9995 

0*0000 

0.0000 

0*0000 

76*9959 

88.7441 

88.7461 

0.0020 

29*9995 

29.9995 

0.0000 

77*9939 

78*9959 

88.1474 

87.3795 

88*1494 

87.5816 

0*0020 

0.0020 

29*9995 

29*9995 

29.9993 

29*9995 

0*0000 

0.0000 

79*9959 

80*9959 

81.9959 

87*0301 

86.5674 

86.1389 

87*0522 

86*5695 

86.1403 

0-0020 

0*0020 

0-0013 

29*9995 

32*9994 

32*9994 

29*9995 

32*9994 

32-9994 

0*0000 

0*0000 

0-0000 

82*9939 

85*7709 

85.7723 

0.0013 

32*9994 

32*9994 

0-0000 

83*9939 

84*9959 

65.4701 

83.2401 

85*4708 

85.2413 

_  0*0006 

0-0013 

32*9994 

32*9994 

32*9994 

32*9994 

0*0000 

0*0000 

85*9959 

86*9959 

87*9939 

88*9959 

89*9959 

90*9939 

91.9939 

92*9959 

93.9959 

85 *0856 

85*0087 

85*0101 

85.0904 

85*2484 

85*4811 

85,7846 

86*1534 

86,5839 

8S*0863 

85*0087 

85.0101 

83.0897 

85*2477 

85*4798 

85,7833 

86.1534 

86,3846 

0*0006 

o.oooo 

0.0000 

-0.0006 

-0*0006 

-0*0013 

-0,0013 

32*9994 

32*9994 

35.9994 
35*9994 

35*9994 

35*9994 

35.9994 
35*9994 
35*9994 

32*9994 

32.9994 

35.9994 
35.9994 

35.9994 

35-9994 

35.9994 

0*0000 
*  0*0000 
0*0000 
0*0000 
0*0000 
0*0000 
0.0000 

-0*  0020 
-0.0013 
-0.0020 
-0.0020 
-0*0020 
-0-002? 

=  -0*0020 
-0.002? 

35*9994 

35.9994 

35.9994 

38.9994 
38*9994 

38*9994 

38*9994 

38*9994 

0*0000 

0.0000 

0*0000 

0*0000 

0*0000 

0*0000 

0*0000 

0.0000 

94*9959 

_  95.9939 

96*9939 

97*9959 

98*9959 

99*9959 

87.0707 

87.6022 

88*1707 

88*7688 

89.3854 

90*0123 

87.0687 

87,6001 

88*1687 

88*7660 

89,3833 

90*0096 

35*9994 

38,9994 

38*9994 

38*9994 

38*9994 

38*9994 

2.  Sample  decimal  on-line  output.  (See  Fig.  1,  Sine  and  Step  functions  run 
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Sample  input  deck  to  7094  print  program. 


SINE  AND  STEP  FliNCTIONS 
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Sample  output  of  7094  print  program. 
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Azimuth  angle  (sine  function). 


SINE  RND  STEP  FUNCTIONS 
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Azimuth  error  (sine  function). 
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Elevation  angle  (step  function). 


SINE  RND  STEP  FUNCTIONS 
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Fig.  10.  Elevation  error  (step  function). 


Hinwizb 

COOLW 

OOOGIO* 

000600'  OOOSOO' 

OOOlOO* 

000600 

yoyy3  Hinwizb 


76 


Fig.  11.  Azimuth  angle  (5th  order  polynomial). 
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Fig.  12.  Azimuth  angle  (5th  order  polynomial).  Plotted  with  different  time  limits. 
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Fig.  13.  Elevation  angle  (impulse  function) . 
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Fig.  14.  Elevation  error  (impulse  function). 
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